Preparing topological PEPS on a quantum computer 
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Simulating of exotic phases of matter that are not amenable to classical techniques is one of the 
most important potential applications of quantum information processing. We present an efficient 
algorithm for preparing a large class of topological quantum states - the G-injective Projected 
Entangled Pair States (PEPS) - on a quantum computer. Important examples include the resonant 
valence bond (RVB) states, conjectured to be topological spin liquids. The runtime of the algorithm 
scales polynomially with the condition number of the PEPS projectors, and inverse-polynomially in 
the spectral gap of the PEPS parent Hamiltonian. 



Creating and studying exotic phases of matter is one of 
the most challenging goals in contemporary physics. The 
increasingly sophisticated simulation abilities of systems 
such as cold atoms in optical lattices, trapped ions or 
superconducting qubits make this challenge accessible by 
means of Feynmann's original idea of using highly control- 
lable quantum systems in order to simulate other quantum 
systems. Among those exotic phases, non-abelian topo- 
logically ordered states and topological spin liquids - such 
as resonating valence bond (RVB) states in frustrated 
lattices - are probably the holy grails of this area of quan- 
tum state engineering. Progress on the creation of such 
exotic phases in various experimental systems has accel- 
erated rapidly in recent years, including cold atoms [1], 
ion traps [2], photonic devices [3] and superconducting 
devices [4]. 

Recently [5] , a very general way of constructing quan- 
tum states on a quantum computer was proposed. The 
wide applicability of the method lies in the fact that there 
is a variational class of quantum states, called Projective 
Entangled Pair States (PEPS), which has a simple local 
description but is nonetheless complex enough to approx- 
imate the low-energy sector of local Hamiltonians. (A 
review of the analytical and numerical evidence for this 
can be found in [6] and the references therein.) However, 
a crucial technical assumption in the main result of [5], 
called 'injectivity', excludes any possibility of constructing 
quantum states with topological order. 

The main aim of this article is to significantly extend 
the result of Ref. [5] to include exotic topological quantum 
phases, by proving: 

Main result For any finite group G, a G-injective PEPS 
can be prepared on a quantum computer in polynomial 
time, when the inverse ground state gap of the associated 
parent Hamiltonian scales at most polynomially in the 
system size. 

'G-injectivity', introduced only recently in [6] (and 
explained more fully below), is a substantially weaker 
requirement than injectivity, which explicitly allows for 



topological order. A compelling example of the signifi- 
cance of this result is the very recently proven fact that 
the resonating valence bond (RVB) state in the Kagome 
lattice (conjectured to be a topological spin liquid), is a 
Z2-injective PEPS, with numerical evidence that the gap 
assumption is also verified [7]. Our result therefore gives 
one way in which the RVB state (and other topological 
states) can be prepared efficiently on a general quantum 
simulator. Engineering exotic quantum states by quan- 
tum simulation complements research aimed at finding 
materials that directly exhibit topological behaviour, and 
is already beginning to bear fruit experimentally [1—4]. 

In the following section, we summarise basic notions 
of PEPS required in this work, and introduce the class 
of G-injective PEPS which includes many of the impor- 
tant topological quantum states. We then briefly review 
the algorithm of Ref. [5] for preparing injective (non- 
topological) PEPS, before proceeding to show how this 
algorithm can be extended to the much larger class of 
G-injective PEPS, thereby allowing efficient preparation 
of many exotic topological quantum states. Finally, we 
close with some concluding remarks and open questions. 

Projected Entangled Pair States For simplicity, we 
will focus in this letter on PEPS defined on a square lat- 
tice, but the results can be generalized to other lattices. 
An (unnormalized) PEPS can be described as follows. 
Place maximally entangled states of dimension D along 
all edges of the lattice. To each vertex v, apply a lin- 
ear map A v : ((D^)® 4 -> <C d in the four £>-dimensional 
systems labeled by I, t, r, b (for 'left', 'top', 'right' and 'bot- 
tom'), where A" = Y,i-i t,r b ^-i-Urb H)(ltrb\, The resulting 
vector in (<C d )® N is the unnormalized PEPS, N being the 
number of vertices in the lattice. For the purposes of this 
work, since local unitaries do not change the complexity 
of preparing a state, by taking the polar decomposition 
of A we can assume without loss of generality that A is 
positive-semidefinite. When A is invertible, we call the 
PEPS injective [6]. 

A particularly interesting class of PEPS is the class 
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of G-isometric PEPS, defined for any finite group G as 
follows. Take a semi-regular representation of G [6] - that 
is, a representation U g = ® a V^ ® l ra having at least one 
copy of each irrep a. Note that the regular representation 
is exactly the one for which r a is the dimension d a of the 
irrep Vg for all a. We can define the re-weighting map 

A = ® Q (^) 4 l da ®l ra (i) 

which is real, diagonal, commutes with U g and satisfies 
Tr A 4 U g = \G\5g te . (For the regular representation A = 1.) 
The PEPS is then defined by taking, for all v. 

A " = JWl E A[7 9 ® AU 9 ® A Ug 9 AU g . (2) 
1^1 geG 

G- isometric PEPS were originally defined in [6] only for 
the regular representation, and shown in that case to be 
exactly the quantum double models of Kitaev [8]. Here, we 
generalise the definition of G-isometric to any semi-regular 
representation [9] . If, on top of a G-isometric PEPS, we ap- 
ply a further invertible (and w.l.o.g. positive-definite) lin- 
ear map A" : <E d -> <D d , we obtain a 'G-injective' PEPS [6]. 
(Here, d is the dimension of the symmetric subspace asso- 
ciated with the group.) The parallel with plain injective 
PEPS is clear. Both are defined by invertible maps on top 
of a G-isometric PEPS. In the case of injective PEPS, the 
group is the trivial one and the representation is simply 
Id (d copies of the left-regular representation of the trivial 
group). 

G-isometric PEPS have very nice properties, coming 
from their topological character, which are inherited by 
the more general G-injective PEPS. For instance, for each 
G-isometric PEPS \ip) there exists a local frustration- 
free Hamiltonian (called the PEPS "parent Hamilto- 
nian" [(>]), consisting of commuting projectors and hav- 
ing as ground space the subspace (over-) spanned by 
{\ip;K):K=(g,h),[g,h]=0}. (Here, \ip; K) is the PEPS 
obtained by the same maps A, except that we first apply 
an additional U® v to exactly one vertical strip V and 
U® H to exactly one horizontal strip H in the initial collec- 
tion of maximally entangled states [6]). This generalises 
to G-injective PEPS, except that the local Hamiltonian 
terms are no longer necessarily commuting projectors. 

We will denote by l^ 1 -^*) the G-injective PEPS de- 
fined by applying the map A 3 to vertex j for j = 1, . . . , t 
(and identity to the rest of the vertices) on top of the 
G-isometric PEPS, and define the states |A — .A*; K) anal- 
ogously to above, which again (over-) span the ground 
space of a frustration- free local parent Hamiltonian H t . 

Preparing injective PEPS We first briefly review 
the algorithm of [5] for preparing injective PEPS on a 
quantum computer. Let H t be the parent Hamiltonian of 
the partially constructed state ^•••A*}. The algorithm 
starts at t = with maximally entangled states between 
all pairs of adjacent sites in the lattice, and proceeds by 



successively projecting onto the ground states of H t for 
t = l...N until the final state \A 1 ---A N ) is reached. 

Since the ground state P t of H t is a complex, many- 
body quantum state, it is not immediately clear (i) how to 
efficiently perform the projective measurement {P t ,Pj-} 
onto the ground state. Furthermore, measurement in 
quantum mechanics is probabilistic, so even if this mea- 
surement can be performed, it is not at all clear (ii) how 
to guarantee the desired outcome P t . 

The answer to (i) is to run the coherent quantum phase 
estimation algorithm [10, If] for the unitary generated by 
time-evolution under H t . (Time-evolution under the local 
Hamiltonian H t can be simulated efficiently by standard 
Hamiltonian simulation techniques [12].) If £/. a k \tpk} is 
the initial state expanded in the eigenbasis of H t , then the 
phase estimation entangles this register with an output 
register containing an estimate of the corresponding eigen- 
value: Y,k a k IV'fe) \Ek)- Performing a partial measurement 
on the output register to determine if its value is less than 
At (the spectral gap of H t ) completes the implementation 
of the measurement {P t ,P^}. (See [5] for full details.) 

The solution to (ii) is more subtle, and makes use of 
Camille Jordan's lemma of 1875 on the simultaneous block 
diagonalization of two projectors, which we first recall 
(see also the related CS decomposition [13]): 

Lemma 1 (Jordan [14]) Let R andQ be two projectors 
with rank s r = rank R and s q = rank Q respectively. Then 
both projectors can be decomposed simultaneously in the 
form 

R=®R k g = ©Q fc , (3) 

fe=l fe=l 

where Rk,Qk denote rank-1 projectors acting on one- or 
two-dimensional subspaces. The eigenvectors |rfc) ,|r^) 
and |qfc) , |<7jt) °f the 2x2 projectors Rk and Qk are related 
by 

Vk) =\fdk\q k ) + y/l-d k \qi) 
H) = -Vl :r d~k\qk) + Vd~k\q 1 k) 

k) =Vd~k\r k }-^l-dkH) 
l?fe)= \f^ r d~kVk) + \fd~k\ri) . 

Ref. [5] shows that if the current state is in the block con- 
taining the ground state of H t , then the PEPS structure 
guarantees that the probability of a successful projection 
onto P t+ i is lower-bounded by K,(A t+1 )~ 2 , where n(A t+1 ) 
is the condition number of the matrix A t+1 . Assume for 
induction that we have already successfully prepared the 
(unique) ground state of H t . We first attempt to project 
from this state onto the unique ground state of H t +\ by 
measuring {P M ,P^ +1 }. If this fails, we attempt to project 
back to the state we started from by measuring {P t ,P^-}, a 
technique introduced by Marriott and Watrous [15] in the 
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context of QMA-amplification. If this "rewind" measure- 
ment succeeds, then we're back to where we started and 
can try again. What if the "rewind" measurement fails? 
By Lemma 1, we can only be in the excited state from 
the same block, so we can still try to project "forwards" 
with the same lower bound on the success probability. 
Thus iterating forwards and backwards measurements 
until success generates a Markov chain with successful 
projection onto the ground state of H t+ \ as the unique 
absorbing state. Moreover, since the success probability 
in each step is bounded away from zero, this converges 
rapidly to the desired state, allowing us to move from 
the ground state of H t to the ground state of H t +i in 
polynomial time. 

Preparing G-injective PEPS Consider the algo- 
rithm of the preceding section from the perspective of 
G-injective PEPS. An injective PEPS can always be 
viewed as a G-injective PEPS for the representation 1 
of the trivial group. The algorithm starts from the state 
consisting of maximally-entangled pairs between each site, 
and transforms this into the desired state by projecting 
onto the ground states of a sequence of injective parent 
Hamiltonians. But the initial state is none other than the 
G-isometric PEPS corresponding to the representation 1 
of the trivial group. This hints at a generalisation of the al- 
gorithm to G-injective PEPS for arbitrary groups G: start 
by preparing the corresponding G-isometric PEPS, and 
successively transform this into the desired G-injective 
PEPS by projecting onto the ground states of the sequence 
of G-injective parent Hamiltonians (see Table I). 

Input: G-injective A" denned on an iV-vertex lattice; e > 0. 
Output: k/>) e spanIA 1 , . . . , A l ;K) with probability > 1 - e. 

1. Prepare corresponding G-isometric PEPS 

2. For t = 1 to N 

2.1. Measure {Pt+i,Pt+i} on \ip). 

2.2. While measurement outcome is P£ +1 

2.2.1. Measure {P t ,P^}. 

2.2.2. Measure {P M , Pt +1 } . 

TABLE I. Preparing a G-injective PEPS. Ht (t = 1 - . . N) is 
the parent Hamiltonian for the G-injective PEPS \A* . . . A*), 
Pt the projector onto its ground state subspace. Note that 
by specifying A v , we are implicitly selecting a particular semi- 
regular representation of G. 

There are, however, two obstacles to implementing 
this approach, (i) The initial G-isometric PEPS can be 
a substantially more complicated many-body quantum 
state than the trivial product of maximally-entangled 
pairs we must prepare in the injective case, (ii) Since 
G-injective parent Hamiltonians are topological, they have 
degenerate ground state subspaces. But the Marriott- 
Watrous "rewinding trick" [15] relies on the measurement 
projectors being 1-dimensional; it breaks down in general 



for higher-dimensional projectors. 

It turns out that there is a direct solution to (i) . Ref. [6] 
proves that, for any group G, the parent Hamiltonian 
of the G-isometric PEPS for the regular representation 
corresponds precisely to a quantum-double model [8, 16]. 
But Ref. [17] shows that ground states of quantum-double 
models can be generated exactly by a polynomial-size 
quantum circuit. We can therefore use this circuit to 
efficiently prepare the G-isometric PEPS for the regular 
representation of G. What if the representation we require 
is semi-regular? In fact the G-isometric PEPS for any 
semi-regular representation is equivalent to the one for 
the regular representation (up to local isometries), by 
simply regrouping the tensors. (See the appendix for a 
proof.) 

The second obstacle is more delicate. As described 
above, the Marriott-Watrous "rewinding trick" used in the 
injective case [5] works because, thanks to the Hamiltoni- 
ans H t and H t +\ at each step having unique ground states, 
there is only ever one 2x2 block involved in the back-and- 
forth measurements. However, in the G-injective case, 
the Hamiltonians no longer have unique ground states, 
and there are multiple 2x2 blocks corresponding to dif- 
ferent ground states. Thus when we "rewind" a failed 
measurement, the backwards measurement could project 
us back into any superposition of states from any of the 
blocks corresponding to the ground state subspace. Now, 
A t+l is only invertible on the G-symmetric subspace, so it 
necessarily has some zero eigenvalues. Hence k{A 1+1 ) = oo 
and the lower bound K,(A t+1 )~ 2 = on the probability of 
a successful forward measurement is useless. Although 
there may still exist some ground state \ipl) of H t which 
has positive probability of successful forward transition 
to a ground state of H t+ i, this does not rule out existence 
of another ground state |*0f) of H t for which the proba- 
bility of a successful forward transition is 0. In the worst 
case, if a forward measurement fails and we end up in a 
state |vt + i), the rewinding step could have probability 1 
of transitioning back to \4>t), so that we remain stuck 
forever bouncing back and forth between \ip 2 ) and |<£>t + i)- 

To overcome this, we must show that if we start from 
the G-isometric state, then the structure of G-injective 
PEPS ensures that this situation can never occur. To 
prove this, we need the following key lemma: 

Lemma 2 Let P t and Pt+i denote two projectors on 
the ground state subspace of the partial PEPS par- 
ent Hamiltonains H t and H t +\ for \A 1 . . . A 1 ) and 
\A X . . .A\ A t+1 ). The overlap d k between P t and P t+1 (cf. 
Lemma 1) is lower-bounded by d m i D > K(A t+1 \s G )~ 2 , where 
K(A t+1 \s G ) ■= a ma , x (A t+ i\s G )/o- min (A t+ i\s G ) is the condi- 
tion number restricted to the G-symmetric subspace Sg- 

Proof The minimum overlap d m i n between projectors Pt 
and Pt+i is given by 

d m m = min max | ( ip t I ipt+ 1 ) 1 2 , (4) 
m) m+i) 
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where \ipt) and \ipt+i) are states in the respective ground 
state subspaces kerif t and kerH t +i. 

Now, ker H t is spanned by the partially constructed 
PEPS IA 1 ,..., A 1 ; K), with different boundary conditions 
K giving different ground states). Thus we can decompose 
any \ip t ) s ker77 t as a linear combination of partial PEPS: 
\iPt) = Zc k \A\...,A t ;K k ). 

Ht+i is obtained from H t by replacing all the 
G-isometric local Hamiltonian terms at one vertex 
with the G-injective terms. So applying A t+1 to any 
\A , . . . , A l \K) takes us to the next ground state subspace. 
Therefore, the state 



IVt+i) 



A t+ i\i> t ) 



yJ(1>t\Al +1 A t+1 \tl> t ) 



(5) 



is contained in kerff t+ i. Choosing \ipt+i) = l¥>t+i) m 
Eq. (4), we obtain the lower bound 

A -s • 1/ / I \l 2 -s • \^t\M+\ \^t)f 

etmin > m.m\(ip t \tp t+1 )\ > mm r . (6) 



I'M 



1'W ^ t \A\ +1 A t+1 \ij t ) 



It is immediate from the definition of G-injective PEPS 
that the ground states of H t are symmetric, so that the 
projector P t is supported on the symmetric subspace 
Sq. Thus the minimisation is over symmetric states 
and, recalling that w.l.o.g. A t is positive-semidefinite, we 
obtain the claimed bound 

. {^ t \ A t+1 \ s \^ t f a min (A t+1 \ s ) 2 

dmin > mm 2 r—r- > ^ , (7 

by the variational characterisation of eigenvalues. □ 

Runtime We are now in a position to establish the 
runtime of the algorithm given in Table I. We start by 
bounding the failure probability of growing the partial 
PEPS by a single site. 

Lemma 3 The measurement sequence depicted in Fig. 1 
with the two projective measurements {Pt,P^} and 
{Pt+i, Pt + ±} has a failure probability bounded by 



Pfaii(m) < 



2d„ 



(8) 



after m-subsequent measurement steps, where c? m i n = 
mhifc d k is the minimal overlap between the eigenstates of 
Pt and Pt+\- 

Proof Let Q x = P f+1 , Q = P^ and fl x = P t , Rq = P t \ 
in accordance with the notation in Lemma 1. Hence 
Qi projects on to the new ground state subspace, 
whereas the i?i is the projector on to the old ground 
state subspace. If we start in some state \ip) = R\ \ip), 
the probability of failure of the measurement se- 
quence depicted in Fig. 1 after m steps is Pf a n(m) = 




FIG. 1. The sequence of outcomes of the binary measurements 
{Pt,Pt} and {Pt+i,Pi+i}- We are initially in an eigenstate 
\ipt) of the projector P t , and want to transition to a state 
in the subspace Pt+i- With non-zero probability, the first 
{Pt+i, Pt+i} succeeds with outcome Pt+i- If it fails, we have 
prepared a state in the P£ +1 subspace. We "unwind" the mea- 
surement by measuring {Pt,Pt} again. Upon repeating the 
{Pt+i, Pt+i} measurement, we again have a non-zero probabil- 
ity of successfully obtaining the P t+1 outcome. If we fail again, 
we repeat the procedure until success. 



S sl ,..., Sm Tr(Q R Sm Qo ■ ■ ■ R Sl Qo \i/>){1>\ Q R S1 ■ ■ ■ QoRsM- 

Note that [QoR s Qo, QoRpQo] = for all s,p. We can 
therefore rearrange this to express pi a i\(m) as the sum 

E (7) M (QoRoQv) 2k (Q RM 2{m - k) W 

k=Q ' fc ' 



((QoRoQo) 2 + (QoRiQo) 2 ) 1 



If we work in the eigenbasis of Qx, the individual 2x2 
block matrices take the form 



Q " = (o o) ' ^ = (v^fr 



Vdfc(l-dfc)^ 



l-d k 



(9) 



Since \^} is left invariant by R\, we have that 
IV') = T.k c k\i~k), where in this basis every \r k ) = 
(\/d~k \/l - dk) T by Lemma 1. We are therefore left 
with 

Pfaii(m) = £|c fc | 2 (1 - d k ) (1 - 2d fe (l - d k )) m , (10) 
fc 

with d k e [0,1] and Z fc |c fe | 2 = 1. 

Since (1 -x) < e~ x , we may bound (1 - d k )(l - 2d k (l - 
d k )) m < (1 - d k )e- 2mdk ^- dk \ Furthermore, we have that 
(l-d k )e- 2mdk{1 - dk) < l/2md k by Taylor expansion. If we 
now choose the largest factor (2mdfc) _1 < (2md m i n ) _1 , we 
can bound the total failure probability by Eq. (8). □ 

We use this to bound the overall runtime. 

Theorem 4 (Runtime) Let A" be G-symmetric tensors 
defining a PEPS on an IS! -vertex lattice. A state in 
the subspace spanned by the corresponding G-injective 
PEPS \ A 1 . . . A N ; K) can be prepared on a quantum com- 
puter with probability 1 - e in time 0(N 4 K 2 jA^ 1 e^ 1 ) , 
with additional classical processing 0(Nd 6 ), where A = 
min t (A t ) is the minimal spectral gap of the family of 
parent Hamiltonians H t for \A l . . . A f ) (t = 1 . . . N), and 
k g = max t K(A*|s G ). 
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Proof The algorithm in Table I first prepares the initial 
G-isometric PEPS, which can be done exactly in time 
0(N log N) [17], and then transforms this step by step 
into the G-injective PEPS, with one step for each of the 
N vertices of Q. Each step has a probability of failure 
Ptailijn) if we repeat the back- and- forth measurement 
scheme m times. We need to ensure that the total success 
probability is lower-bounded by (1 - Pia,i\{m)) N > 1 - e. 
Since (1 - x) N > 1 - Nx, we can use Lemma 3 to bound 



(1 - Rail (m)) N > 1 



N 



2md„ 



(11) 



so we want N/2md m i n < e. Since d m i n > Kg 2 by Lemma 2, 
we choose m > Nkq/2€ at each step. We therefore need 
to perform 0(N 2 KqE^ 1 ) quantum phase estimation pro- 
cedures, each of which has runtime 0(N 2 / A -1 ) to ensure 
that we are able to resolve the energy gap of the parent 
Hamiltonian [ ]. (Note that the notation O(-) suppresses 
more slowly growing terms such as exp(y / ln N/ A).) The 
classical bookkeeping required to keep track of the Hamil- 
tonians is the same as in [ ]. Putting all this together, we 
arrive at the total runtime stated in the theorem. □ 

Discussion We have shown how the Marriott- Watrous 
rewinding technique combined with the unique structure 
of G-injective PEPS can be used to transition from one 
state to the next (see Table I), successively building up 
the desired quantum state even when that state has topo- 
logical order and the ground states are degenerate. There 
are a number of alternative techniques that could poten- 
tially be used to achieve the same thing. In each case, 
the key to proving an efficient runtime is our Lemma 2. 
In many cases the existing results in the literature as- 
sume non-degenerate ground states, so would need to be 
generalised before they would apply to the topologically 
degenerate ground states considered here. 

Standard adiabatic state preparation has the disadvan- 
tage that it requires a polynomial energy gap along a 
continuous path joining the initial Hamiltonian with the 
final one. But the "jagged adiabatic lemma" of Ref. [19] 
shows that such a path connecting a discrete set of gapped 
Hamiltonians always exists if the ground states are unique 
and each ground state has sufficient overlap with the next. 
The latter is precisely what we prove in Lemma 2. For 
the 'injective' case of [5], this is sufficient to show that 
adiabatic state preparation is an efficient alternative to 
the "rewinding trick" . Our results suggest that the jagged 
adiabatic lemma could be generalised to the case of non- 
unique ground states. 

More general are the methods of [2( )] , which subsume 
the jagged adiabatic lemma and the Marriott- Watrous 
technique. The results in [20] do not immediately apply to 
degenerate ground states, but if they can be generalised 
they could potentially improve the polynomial depen- 
dence on the required error probability to a logarithmic 
one. Another potential alternative is the recent quantum 



rejection sampling technique of [21], which gives quadratic 
improvement over Marriott- Watrous rewinding by a clever 
use of amplitude amplification. Finally, the spectral gap 
amplification technique of [22], which cites injective PEPS 
preparation [5] as a potential application, may also be 
applicable. In all cases, the techniques would first need 
to be generalised to handle degenerate ground states. If 
this can be done, our Lemma 2 would imply efficiency of 
the resulting algorithm. 

The conditions required for efficient preparation in 
Theorem 4 (inverse-polynomial scaling with system size 
of the spectral gaps of the partial parent Hamiltonians 
and polynomial scaling of the condition numbers of the 
PEPS projectors) are very reminiscent of the conditions 
(local gap and local topological quantum order) required 
for stability of the spectral gap of local Hamiltonians [23] . 
It is also conjectured that the spectral gap of the parent 
Hamiltonian should be closely related to the condition 
number of the PEPS projectors. It would be interesting to 
understand better the relationships between these various 
conditions. 

The technique we introduced, of constructing a com- 
plex many-body quantum state by starting from an easily- 
constructable state and successively transforming it into 
the desired state, is very general. Although we have ap- 
plied it here to G-injective PEPS, as a class of states 
including many important topological quantum states 
such as the RVB state, our algorithm can be generalised 
to other classes of tensor network states, such as string- 
net models [24] and models constructed from Hopf alge- 
bras [25]. 

Acknowledgements DPG and TSC would like to 
thank the hospitality of the Centro de Ciencias Pedro 
Pascual and the Petronilla facility in Benasque, where 
part of this work was carried out. TSC is funded by the 
Juan de la Cierva program of the Spanish science min- 
istry. DPG and TSC are supported by Spanish grants 
QUITEMAD, I-MATH, and MTM2008-01366. MS is sup- 
ported by Austrian SFB project FoQuS F4014. KT is 
funded by the Erwin Schrodinger fellowship, Austrian 
Science Fund (FWF): J 3219-N16. FV is supported by 
EU grants QUERG and by the Austrian FWF SFB grants 
FoQuS and ViCoM. 



[1] I. Bloch, J. Dalibard, and S. Nascimbene, Nature Physics 
8, 267 (2012). 

[2] R. Blatt and C. F. Roos, Nature Physics 8, 277 (2012). 
[3] A. Aspuru-Guzik and P. Walther, Nature Physics 8, 285 
(2012). 

[4] A. A. Houck, H. E. Tureci, and J. Koch, Nature Physics 
8, 292 (2012). 

[5] M. Schwarz, K. Temme, and F. Verstraete, Phys. Rev. 

Lett. 101, 110502 (2012). 
[6] N. Schuch, I. Cirac, and D. Perez-Garcfa, Annals of 



Physics 325, 2153 (2010). 

[7] D. Poilblanc, N. Schuch, D. Perez-Garra'a, and J. I. Cirac, 
Physical Review B 86, 014404 (2012). 

[8] A. Kitaev, Annals of Physics 303, 2 (2003). 

[9] The generalisation of 'G-isometric to any semi-regular 
representation is justified by the fact that, for a given G, 
all G-isometric PEPS are equivalent (up to isometries) 
to the G-isometric PEPS for the regular representation 
by decomposing the tensor in a different way. A proof is 
given in the Appendix. 
[10] M. Nielsen and I. Chuang, Quantum Computation and 
Quantum Information (Cambridge University Press, Cam- 
bridge, 2000). 

[11] E. Knill, G. Ortiz, and R. D. Somma, Physical Review 
A 75, 012328 (2007). 

[12] D. Berry, G. Ahokas, R. Cleve, and B. Sanders, Com- 
munications in Mathematical Physics 270, 359 (2007), 
10.1007/s00220-006-0150-x. 

[13] G. H. Golub and C. F. van Loan, Matrix Computations, 
3rd ed. (Johns Hopkins University Press, 1996). 

[14] C. Jordan, Bull. Soc. Math. France 3, 103 (1875). 

[15] C. Marriott and J. Watrous, Computational Complexity 
14, 122 (2005). 

[16] N. Schuch, D. Perez-Garcia, and I. Cirac, Phys. Rev. B 

84, 165139 (2011). 
[17] M. Aguado and G. Vidal, Phys. Rev. Lett. 100, 070404 

(2008). 

[18] A. Harrow, A. Hassidim, and S. Lloyd, Physical review 
letters 103, 150502 (2009). 

[19] D. Aharonov and A. Ta-Shma, SIAM Journal on Com- 
puting 37, 47 (2007). 

[20] S. Boixo, E. Knill, and R. Somma, (2010), 
arXiv:1005.3034. 

[21] M. Ozols, M. Roetteler, and J. Roland, in Proceedings 
of the 3rd Innovations in Theoretical Computer Science 
Conference, ITCS '12 (ACM, New York, NY, USA, 2012) 
pp. 290-308. 

[22] R. Somma and S. Boixo, (2011), arXiv:1110.2494. 
[23] S. Michalakis, (2012), arXiv:1206.6900. 
[24] N. Schuch, (private communication). 
[25] O. Buerschaper, M. Mombelli, M. Christandl, and 
M. Aguado, (2010), arxiv:1007.5283[cond-mat]. 



7 



Appendix: G-isometric PEPS 

In this appendix, we use the argument described in [G] 
for the Toric Code and RVB states, generalised here to 
arbitrary G-isometric PEPS, to see that the G-isometric 
PEPS for any semi-regular representation is equivalent to 
the one for the regular representation. Let us start with 
a semi-regular representation U g of a group G, and let 

B = r^E A ^9® AU g ®AU g . (12) 



indices, is given by 

(13) 

By calling g = g 1 1 g[ and using Eq. (1) it is not difficult 
to see that 

{91929394.1 C ^ C \9i9 25354) 

\^\ r=l r=3 

(14) 



We will show how B can indeed be seen as the G-isometric equals 1 if and only if there exigt g such that ^ = g , for 

PEPS corresponding to the regular representation - pos- all % _ otherwise, the expression is identically zero, 

sibly composed with an isometry which embeds the initial Therefore = (|C|) 1 £ i?f for the regular repre- 

Hilbert space into a sufficiently large one. As explained sen tation R g , hence the new PEPS C is the G-isometric 

in the main text, the latter can be prepared efficiently on pEps corresp onding to the regular representation, 
a quantum computer by other means. 




FIG. 2. (a) illustrates the decomposition of the original tensor 
in the tensors A. We mark in white the bonds in which we have 
U g and in black those in which we have U g . (b) illustrates the 
new way of grouping the tensors to get a G-isometric PEPS, 
called C. The bonds of this new tensor are numbered clockwise 
as in the figure. 



To show this, we decompose the tensor B into two 
tensors of the form A = (y^Wiy 1 Z g AU g <S> AU g ® \g) 
(where U g and U g are interchanged as needed, as shown 
in Fig. 2(a)). By regrouping these new tensors, we obtain 
a new PEPS decomposition of the same state, where now 
the bond dimension is \G\ (Fig. 2(b)). The resulting tensor 
C (Fig. 2(c)), as a map from the virtual to the physical 



